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Abstract 



> 

t^J" ' The central solar temperature T and its uncertainties are calculated in 

helioseismologically-constrained solar models. From the best fit to the con- 
I/-) ■ vective radius, density at the convective radius and seismically determined 

helium abundance the central temperature is found to be T = 1.58 x 10 7 K, in 
| excellent agreement with Standard Solar Models. Conservatively, we estimate 

that the accuracy of this determination is AT/T = 1.4%, better than that in 
Q_i' SSM. Neutrino fluxes are calculated. The lower limit to the boron neutrino 

O ' flux, obtained with maximum reduction factors from all sources of uncertain- 

ties, is 2a higher than the flux measured recently by SuperKamiokande. 



Helioseismology allows us to look into the deep interior of the Sun, probably more effi- 
ciently than neutrinos (for reviews see |1|-|5|]). The highly precise measurements of frequencies 
and the tremendous number of measured lines enable us to extract the values of sound speed 
and density inside the Sun with accuracy better than 1%. Recently it was demonstrated 
that a comparable accuracy can be obtained for the inner core of the Sun (R/R Q < 0.1) 



Helioseismic data are in agreement with recent Standard Solar Model (SSM) calculations, 
which use accurate equations of state, recent opacity tables and include helium and heavier 
elements diffusion [p|,[7|p|dTl[1 , see also Ref. fll^] . These SSMs yield central temperatures 
Tssm which differ from each other by not more than 1%. However the uncertainties in 
input parameters, mainly the opacity k and the heavy elements abundance Z/X, result in 
(AT/T) SSM « 1 - 2%. ' 

From helioseismological observations one cannot determine directly temperature of the 
solar interior, as one cannot determine the temperature of a gas from the knowledge of the 
sound speed unless the chemical composition is known. However, it is possible to obtain the 



1 



range of allowed values of the central temperature T, by selecting those solar models which 
are consistent with seismic data. [] 

In this paper we shall address two related problems: accuracy in the determination of 
the central temperature in helioseismologically-constrained solar models and neutrino fluxes 
in these models. 

Our calculations are not model- independent, but we shall use in principle a wider class 
of models in comparison with SSMs, which we call helioseismologically-constrained solar 
models (HCSM). These models are based on the same equilibrium and evolution equations 
as SSMs, but they differ in the choice of some input parameters. We generate this class of 
models by using the FRANEC code for the SSMs and varying the input parameters. 
Each choice of the set of parameters gives some value of T. 

We obtain the range of allowed values of T by selecting those solar models which are 
consistent with seismic data. More specifically, we shall determine the central temperature 
Thcsm, as that of the model which gives the best fit to the seismic data and the uncertainties, 
AT HC sm, corresponding to the range spanned by models consistent with these data. 

As a method of calculation, we shall use the scaling approximation of the results of solar 
models evaluated with FRANEC. Namely we assume that every physical observable Q (e.g. 
temperature, photospheric helium abundance, radius of convective zone etc) depends on the 
input parameters P as 

_o__/_p_\— t 

QsSM V PsSM ) 

where the subscript SSM is related to the (arbitrarily chosen) reference SSM. In this paper 
we use the new determination of the coefficients cxq^p obtained with the latest version of 
FRANEC JK|, which includes diffusion of all elements. 

The precise value of the temperature is governed by four quantities P: the radiative 
opacity k, the fraction of heavy elements Z/X, the astrophysical factor of the p+p reaction 
Six, and the solar age t & . If these four quantities are rescaled with respect to the values 
used in the SSM calculation by a multiplicative factor p, = Pi/Pi^ssu (Pi=k, z, s and t 
respectively), the central temperature scales as: 

T = T 55M k ai4 z - 078 s- 014 t 084 . (2) 

We shall proceed now to determine the range of the scaling variables, Pi, allowed by 
seismic observations. We remark that the errors on Sxx and age t & are small (1% or less), 
so that they weakly affect T. On the other hand, the uncertainties on k and Z/X are of 
the order of 10% and their influence on T is important. Furthermore, these uncertainties 
do not correspond to clear experimental or observational errors, rather they are determined 
by judicious comparison among published values. As an example, the uncertainty on k can 
only be estimated from the comparison among recent theoretical calculations. 

We shall thus concentrate on the two most important scaling variables, k and z. 



Here and everywhere below T is the central temperature. 
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We want to emphasize that this not a determination of the above-mentioned parameters 
directly from seismic data, since solar models are involved in this evaluation. Rather we 
intend to determine those solar models which are consistent with helioseismology. Let us 
also remark that we are restricting to uniform, although generous, variations of opacity. 

As seismic "observables" we choose three independent physical quantities determined 
most accuratelly by seismic observations (see Ref. J7|), namely the photospheric helium 
abundance, Y p h, the depth of convective envelope, Rb, and the density, pb, at the bottom 
of convective envelope (see the data in Table 1). A fourth seismic "observable", the sound 
speed at the convective radius is traditionally considered, e.g. |J. We have not included it 
in our list since, as shown in Ref. it is not an independent one. 

By numerical experiments, using FRANEC, we found the following dependence on z and 
k (for dependence on the other variables see Table |) 

7? — f? ^-0.046 ^-0.0084 /o \ 

Rb — ttb,SSM z k [6a) 

Pb = p b ,ssMZ 0A7 ^ 5 (3b) 

Y ph = Y ph , SSM z°- 3l k 0M (3c) 



We remark that p& and Rb have almost equal restriction power. In fact the error in Rb 
is much smaller than that in pb (see Table |) , but the dependence on the scaling parameters 
is much stronger for pb, see Eq. (pT3|) . 

We also remark that T is mainly determined by Y p h- One sees from Eqs. (||]) and ( jScf ) 
that to a good approximation the dependence of T and Y ph on z and k is just through 
rj = z ■ k 2 , so that one can express T as a function of Y ph : 

Thcsm = Tssm {Y p h,Q/Y phtSSM ) ' 2 , (4) 

Let us determine now the values of z and k which give the best fit to all three quantities 
Qo,i, i.e. Y p h, Rb and Rb. For this we minimize the function 

The corresponding value of the central temperature T(z, k) gives our best estimate Thcsm- 
By starting with BP95 model (the "best model with metal and helium diffusion" of Ref. 
H) with a central temperature T S sm = 1.584 x 10 7 K, one arrives at Thcsm — 1-587 x 10 7 K, 
slightly higher but fully consistent with the former value within the uncertainties of SSMs. 
In fact, this shows that BP95 is in good agreement with helioseismic data. We repeated the 
same procedure starting from a few different SSM calculations, on which enough information 
is available to us: FR97 is the "best" model with He and heavier elements diffusion of 
Ref. [10[|; FR96 is another variant of the same model with the Livermore opacities tables 



calculated for just 12 elements 0; JCD is the "model S" of Ref. Q. 



We remark that BP95 and FR96 adopt the same 12 elements composition [131, whereas 



FR97 uses the newest Livermore opacities calculated for 19 elements ||14||. The values found 



are all within the range Thcsm — (1-573 — 1.587) x 10 7 K, slightly tighter than that of SSM 
predictions, T SSM = (1.567 - 1.584) x 10 7 K, see also Fig. 1. 
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By averaging over all available SSM as starting models we obtain the best estimated 
seismic temperature 



T HCSM = 1.58xl0 7 K. (6) 

The range of acceptable temperatures is determined by those values of z and k such 
that each Q(z,k) is in agreement with helioseismology within the estimated uncertainty. By 
using BP95 as a reference SSM, the allowed domain of z and k is shown in Fig. 2, where we 
also show the limiting temperatures. The resulting uncertainty is ±1%. 

A conservative estimate of (AT/T) H csm is obtained by adding the spread due to starting 
with different SSMs. The result is shown together with (AT/T)ssm estimated with the same 
procedure, assuming a 10% uncertainty in k and Z/X: 

(AT/T) HCS m = ±1-4%, (AT/T)ssm = ±2.7% (7) 

We remark that Eq. (|?|) for both uncertainties is derived using a similar and conservative 
approach, of adding up linearly all error sources. 

Should one add errors in quadrature (so that e.g. AY p h/Y p h = 1.4%, see J7j), it gives 

{AT/T) HOSM = ±0.5%, (AT/T) SSM = ±1.7% (8) 

Fig. 1 summarizes what we have found so far: the central temperatures in HCSMs agree 
very well with ones in SSMs, though the uncertainties in the former models are smaller. 

Finally, we shall discuss a question which a concerned reader certainly must ask: why 
do we not use the sound speed profile as a constraining condition? Indeed, in the region 
0.2 < R/R G < 0.6 the accuracy on the isothermal sound speed squared U = V/p, where 
V is pressure, is better than 0.5%. Why this tremendous accuracy, taken as a function of 
distance, does not give the strongest constraints? 

Let us look to this problem quantitatively. At each fractional distance x = R/R & we 
parametrize U as: 

U(x) = U SS m{x) ^ {x) k° k(:c) s as{x) (9) 

The calculated coefficients, a p (x) = d\nU/d\iap, are plotted in Fig. 3. For each parameter 
p the excluded value of Ap is given by 

ApS Z7v < 10 > 

In case of Z/X, for example, the maximum a z is 0.05 (see Fig. 3) and for AU/U ~ 0.005 
we can exclude variations in Z/X of order of 10% or more, which does not improve our 
knowledge. 

As noted in Ref. |I5|, a change in n(p,T) by a multiplicative constant factor can largely 



be compensated by change in the composition, the sound speed profile remaining approxi- 
mately the same in the intermediate region. Actually the opacity coefficients are an order 
of magnitude smaller than the others, so that a 10% uniform variation of k, which affects T 
to the 1.5% level [see Eq. (|2])], cannot be excluded by studying U(x). 
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Let us come over to neutrino fluxes, (z=pp, Be, B). Their dependence on the central 
temperature T is parametrized as: 



^ = ^ssm[^ • (11) 

V 1 SSM / 

From numerical experiments with FRANEC, we found: [3 pp = —0.92, /3 Be = 7.9, 0b = 18. 

We have determined neutrino fluxes by starting with different SSMs, renormalizing their 
predictions to the same temperature Thcsm = 1-58 • 10 7 K and to the same (updated) 
nuclear cross sections. The resulting fluxes and signals, all very close to each other, have 
been averaged to determine the HCSM predictions shown in Table U where the uncertainties 
correponding to (AT/T)hcsm = ±1-4% are also indicated. 

In order to discuss nuclear physics uncertainties, we shortly review the present status for 
three most important cross-sections. The numbers quoted below refer to 3a errors. 

a) 3 He+ 3 He — > 4 He + 2p. This cross-section was measured recently in LUNA experiment 
P~5| ] at energy corresponding to the solar Gamow peak (~ 20 KeV). The uncertainty in 
astrophysical factor is AS^/S^ = ±6%. 

b) 3 He+ 4 He -> 7 Be + 7. The uncertainty is AS34/S34 = ±12%. 

c) p+ 7 Be — > 8 B + 7. S17 is determined from direct measurements with an error of about 
±30% |L7-19|. The indirect measurement [^TJ] gives S17 consistent with the lower limit of 



direct measurements, although the accuracy of this method is unclear pl|-p3 . 

The errors in fluxes due to cross-sections are calculated using the relations between fluxes 
and cross-sections |24],|2"5|j , <3>£ e oc S34 S 3 ^ 2 and $b oc S34 S 3 ^ 2 S17. We remark that the 



(3a) uncertainty on $£ e due to that on S34 is slightly larger than that corresponding to 
temperature. The 30% uncertainty on due to that of S17 exceeds that due to the solar 
temperature. 

One concludes, see Table |T|, that (3a) nuclear physics uncertainties are at least as im- 
portant as the (generously estimated 1.4%) temperature uncertainty. 
In the end we shall discuss the problem of boron neutrinos. 

The HCSM boron neutrino flux (see Table II) is 18a higher than than the combined 
Kamiokande |26] and SuperKamiokande |27| flux, = (2.58 ±0.19) • 10 6 cm" 2 s _1 ). How- 



ever, the uncertainties in the predicted flux are of great importance for this comparison. 

During the last several years an idea of reconciling the predicted boron neutrino flux with 
that measured by Kamiokande was widely discussed |f2"8"|-|3"l"|j . This discussion was inspired 
by the low value of Sn measured in indirect experiments [^0J and by a possible decrease of 
the central temperature T due to collective plasma effects (through opacity [p2p and due to 
abundance of heavy elements Z. These effects, when correlated, can result in the agreement 
between predicted and measured B-neutrino flux. 

Since in our calculations the opacity and heavy elements abundance are constrained by 
seismic observations, the status of this problem has changed. Indeed, diminishing Thcsm 
by 1.4% and S17 by 30% one obtains from the boron- neutrino flux in HCSM the mini- 
mum flux 3.24 ■ 10 6 cm~ 2 s _1 , which is 3.5a higher than the combined Kamiokande and 
SuperKamiokande flux. 

If one accepts a more reasonable 1% reduction of Thcsm, but takes other errors in a 
correlated way, all diminishing the B-neutrino flux (namely S17 by 30% smaller, S34 by 12% 
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smaller and S34 by 6% higher) we obtain = 2.97 • 10 6 cm^ 2 s _1 , i.e. 2.1a higher than 
SuperKamiokande flux. Only in case of largest possible correlated errors the HCSM flux can 
be reconciled with the SuperKamiokande data. 

Therefore, we have now three solar-neutrino problems: (i) the controversy of Homestake 
and SuperKamiokande results, (ii) the Beryllium neutrino problem, as controversy of gallium 
and SuperKamiokande (or Homestake) experiments, and (iii) Boron neutrino problem, as it 
is described above. 

In conclusion, the helioseismologically constrained solar models (HCSM) give the central 
solar temperature in excellent agreement with SSMs, and with smaller uncertainties. The 
boron neutrino problem exists in this class of models. 
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TABLES 



TABLE I. Seismically determined quantities, Q, and the exponents «q,p of the scaling ap- 
proximation given by Eq.(l). 

0© <y-Q,k OQ^z OQs OLQ£ 

Y ph 0.249 (1 ±4.2%) ~~ 0.61 0.31 0.14 0.20 

R b /R & 0.711 (1 ±0.4%) -0.0084 -0.046 -0.058 -0.080 

p b [gr/cm 3 ] 0.192 (1 ±3.7%) 0.095 0.472 0.86 0.85 



TABLE II. Predictions for neutrino fluxes and signals in the CI and Ga detectors from SSM 
and HCSM. Uncertainties corresponding to (AT ' /T)hcsm = ±1-4% are shown (first error) together 
with those from nuclear cross sections (second error). 





BP95 


SSM 

FR97 


FR96 


JCD 


HCSM 




[10 9 /cm 2 /s] 


5.15 


4.49 


4.58 


4.94 


4.81±0.53±0.59 




[10 6 /cm 2 /s] 


6.62 


5.16 


5.28 


5.87 


5.96±1.49 ± 1.93 


CI 


[SNU] 


9.3 


7.3 


7.5 


8.2 


8.4±1.9±2.2 


Ga 


[SNU] 


137 


128 


129 


132 


133±11 ±8 
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FIGURES 



FIG. 1. For a few recent SSM calculations we present : a) the predicted central solar tem- 
peratures Tssm (diamonds) with the conservative uncertainty (thin bars); b) the values Thcsm 
derived by best fit of helioseismic observables at the convective radius (circles) with the uncertain- 
ties (bars) calculated in the same way as for SSM's. Labels indicating solar models are defined in 
the text. 

FIG. 2. In the (z,k) plane we show the position of the best-fit HCSM (diamond) and of the 
models consistent with helioseismic data (crosses) obtained by using BP95 as a starting model. 
The curves labelled by pj,, Y p h and Rb show the helioseismic constraints due to these quantities. 
Curves labelled by the numbers are isotemperature curves in kz plane. Numbers are values of 
T/Tssm- 

FIG. 3. Logaritmic derivates of the isothermal sound speed squared U with respect to Su 
(dot-dashed line), opacity (dashed line) and Z/X (full line). Note that the opacity coefficient have 
been multiplied by a factor 10. 
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